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0. Introduction 

This paper is a continuation of [Z.l]. There, we gave an effective method for putting the Laplacian A of 
a Riemannian manifold (M, g) into a quantum Birkhoff normal form around a non-degenerate elliptic closed 
geodesic 7, and applied it to the calculation and characterization of the wave invariants a 7 fc at 7. The wave 
invariants, we recall, are the coefficients in the singularity expansion 

Trt/(i)=e (i)+E 7 e 7 (i) 

e 7 (i) ~ a 7 _i(t — Lj + iO)^ 1 + J2T=o a 7 fe (* — ^7 + i0) k log(t — L 7 + iO) 

of the trace of the wave group U(t) — e lt ^ at lengths t = L 7 of closed geodesies. The first purpose of this 
article is to extend the methods and results of [Z.l] to general non-degenerate closed geodesies, i.e. to 7 
whose Poincare map P 7 is any symplectic sum of non-degenerate elliptic, hyperbolic, or loxodromic parts. 
Our second purpose is to generalize to the full non-degenerate case the inverse result of Guillemin that the 
quantum normal form coefficients at non-degenerate elliptic closed geodesies are spectral invariants [G.1,2]. 
It will follow that, for metrics with simple length spectra, the classical Birkhoff normal form of the metric 
around any non-degenerate closed geodesic is a spectral invariant of the Laplacian. 

Let us state the results precisely. The first is that the wave invariants a 7 fe in the general non-degenerate 
case are essentially analytic continuations of the expressions obtained in the elliptic case. They may be 
written in the form 

(0.1) a 1 k = T k ,-i{D)-Ch{x)\ x=Pl 

where 

Chtx) = ; % " 

is the character of the metaplectic representation (with a a certain Maslov index) and where Tk,-i{D) is 
an invariant partial differential operator on the metaplectic group Mp(2n, R) which is canonically fashioned 
from the germ of the metric g at 7. The exact expression for Tu.~i{D) will be given in §5 and leads to the 
following characterization of the wave invariants (cf. [Z.l, Theorem A]): 

Theorem I Let 7 be a non- degenerate closed geodesic. Then a 7 ^ = j 7 7; fe(s; g)ds where: 

(i) I 7 -k(s,g) is a homogeneous Fermi-Jacobi-Floquet polynomial of weight -k-1 in the data {yij,yij, D@ y g} 
with \0\ < 2k + 4 ; 

(ii) The degree o// 7; /c in the Jacobi field components is at most 6k+6; 
(Hi) At most 2k+l indefinite integrations over 7 occur in I^kri 

(iv) The degree of J 7; fc in the Floquet invariants [3j is at most k+2. 

The relevant terminology and notation will be recollected in §1. From (0.1) and from the formulae for 
J-k.-i(D) and Ch, we give a rather simple proof (§6) of the following inverse result (strictly speaking, proven 
only for non-degenerate elliptic closed geodesies in [G.2, Theorem 1.4]): 

Theorem II Let 7 be a non- degenerate closed geodesic. Then the entire quantum Birkhoff normal form 
around 7 is a spectral invariant; in particular the classical Birhoff normal form is a spectral invariant. 

We thus have: 
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Corollary II. 1 Suppose (M,g) is a compact Riemannian manifold with simple length spectrum Lsp{M 1 g) 
and with all closed geodesies non-degenerate. Then from Spec(M, g) one can recover the quantum (and hence 
classical) Birkhoff normal forms around all closed geodesies. 

The hypotheses of the corollary are of course satisfied by generic Riemannian metrics (cf. [Kl, Lemma 
4.4.3]). The corollary therefore answers affirmatively the third question in [Z.2, p. 692], which asks whether 
isospectral manifolds in this class of metrics are locally Fourier isospectral near corresponding closed geodesies. 

Let us now briefly discuss the main ideas in the proofs, and in particular the novel aspects caused by 
the hyperbolic and loxodromic parts of P 1 . 

As in [G.l,2][Z.l], the wave invariants at a closed geodesic 7 will be expressed as non-commutative 
residues of the wave group and its time derivatives &tt = L 1 . For simplicity we will often abbreviate L 7 by 
L. Then we have: 

a lk = res D k t e u ^\ t=L := Res s=0 TrD k t ^^ S \ t=L . 

Since res is invariant under conjugation by (microlocal) unitary Fourier integral operators, the a 7 fc's may be 
calculated by putting the wave group into a microlocal (quantum Birkhoff) normal form around 7 and by 
by determining the residues of the resulting wave group of the normal form. 

The primary step in the analysis of the wave invariants a 7 fc is therefore to put A into this microlocal 
normal form around 7. In the case of non-degenerate elliptic closed geodesies, we recall, the normal form 
was a polyhomogeneous function in the (microlocally elliptic) element D s with coefficients in the transverse 
(elliptic) harmonic oscillators 7| = \{P Vj + y]) [Z.l, Theorem B]. In the general non-degenerate case, the 
normal form will involve a greater variety of quadratic normal forms or 'action operators:' in addition to the 
elliptic action operator JJ there can also occur the real hyperbolic action operators /j 1 and complex hyperbolic 
(or loxodromic) action operators I] ' Re ,I c j ' Im . These hyperbolic actions cause several complications to the 
arguments in the elliptic case: First, they have continuous spectra, and so the construction of the intertwining 
operator to the normal form has to be modified in several ways (§3,4). Second, the wave group of the normal 
form has continuous spectrum and this alters the calculation of its residues (§5). Third, the presence of real 
parts in the Floquet exponents of P 7 complicates the process of determining the normal form coefficients 
from the wave invariants (§6). 

To get acquainted with these action operators and the normal form algorithm, let us consider the very 
first step of "linearizing" ^/A around 7 and of putting the 'linearization' 

_^ n n 

into quantum quadratic normal form. Here, n = dimM — 1, the coordinates (s,yj) are the (re-scaled) 
Fermi normal coordinates around 7, D s = -S^ and Kij is the curvature operator g{R{d s ,d yi )d s ,d yj ). The 
linearization C of \fK is a quadratic Hamiltonian and is the Weyl quantization of a quadratic classical 
Hamiltonian (see [Ho III] and [Ho] for background on Weyl quantizations and normal forms for quadratic 
Hamiltonians). Hence its symbol may be conjugated into normal form by an element of W G Sp(2n, R). The 
operator C itself may be put into normal form by conjugating with the metaplectic operator /u(W) quantizing 
W. As will be seen in §1 , this linear symplectic map is the Wronskian matrix W whose columns consist of the 
Jacobi eigenfields of P 7 . The normal form of the linearized \/A is therefore given by 1Z = piW)* Cp(W)*^ 1 . 
In the elliptic case [Z.l] P 7 was a direct sum of rotations, and the quantum normal form of C had the 

form 

H e =D a + ^H a , H a = J2 aj i* 

3=1 

where the spectrum tr(P 7 ) = {e ±Mj }. In the general non-degenerate case the normal form will similarly 
depend on the spectral decomposition of P 1 . Recall that, since P 1 is symplectic, its eigenvalues pj come in 
three types: (i) pairs p, p of conjugate eigenvalues of modulus 1; (ii) pairs p, p~ x of inverse real eigenvalues; 
and (iii) 4-tuplets p, p, p^ 1 p^ 1 of complex eigenvalues. We will often write them in the forms: (i) e ±la \ 
(ii)e ±Xj , (iii) e ±Mj ±l! ^ respectively (with ctj, \j, p,j,Vj G K), although a pair of inverse real eigenvalues 
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{— e ±A } could be negative. Here, and throughout, we make the assumption that P 7 is non- degenerate in the 
sense that 

II^C^l, (V Pi e<7(P 7 ), (m 1 ,...,m 2 „)eN 2 "). 

Each type of eigenvalue then determines a different type of quadratic action, both on the classical and 
quantum levels (cf. [Ho, Theorem 3.1],[Ar]): 



Eigenvalue type 


Classical Normal form 


Quantum normal form 


(i) Elliptic type 
{e i±a } 


P = \a(r?+y 2 ) 


r :=\a{Dl + y*) 


(ii) Real hyperbolic type 
{e ±A } 


I h = 2Xyn 


I h := X(yD y + D y y) 


(iii) Complex hyperbolic 
(or loxodromic type) 


ph,Re = 2n( yi rn + y 2 r) 2 ) 

jchjm = v{yim _ mfli) 


I ch ' Re = KViD yi + D yiVl + y 2 D V2 + D V2 y 2 ), 
I ch > Im = v{yiD y2 -y 2 D m ) 



In the case where the Poincare map P 7 has p pairs of complex conjugate eigenvalues of moduls 1, q 
pairs of inverse real eigenvalues and c quadruplets of complex hyperbolic eigenvalues, the linearized \/A will 
have the form: 

, p 

jch,Im-\ 



n = D s + ^\Y j a i e J 



q c 



fch.Re . T c 



j=i i=i 

The full quantum Birkhoff normal form is then given by the analogue of Theorem B of [Z.l]: 

Theorem B Assuming 7 non-degenerate, there exists a microlocally elliptic Fourier integral operator W 
from the conic neighborhood ofR + j in T*(N y ) to the corresponding cone in T+S 1 in T*(S 1 x W 1 ) such that 

wVaw- 1 = a, + ^E ^i; + E + E ,Re + ^- h ' Im }+ 

j=i j=i j=i 

fch,Re fch,Im\ 



+ 



fe fh fh fch.Re fch,Im 

• ' 1 pl 1 l 1 ■ ' ' ' 1 q 1 1 1 1 1 1 1 



+ 



D s 

Pk+ 1 (if,...,i c c h - Im ) 



+ 



+ ... 



where the numerators Pj{I I, . . . , 1^, i±, . . . , I^ h ' Im ) are polynomials of degree j+1 in the variables (If, . . . , I^ h - Ir - 
and where the kth remainder term lies in the space (B^t^0 2 (^ k+2 _j^ 1 ^ j 

Here, O n ^> r is the space of pseudodiffercntial operators of order r whose complete symbols vanish to 
order n at (y, n) — (0, 0). Thus, the remainder terms are 'small' in that they combine in some mixture a low 
pseudodifferential order or a high vanishing order along 7. 

Some remarks now on the contents and organization of this paper. Since it is a continuation of [Z.l], we 
have tried to avoid duplicating arguments and calculations which are essentially unchanged from the elliptic 
case. Many of the arguments which remain are still quite analogous to the elliptic case and inevitably 
produce a sense of deja-vu. Our excuse for drawing them out to their present length is that it is not apriori 
clear that the arguments of the elliptic case generalize so neatly to the hyperbolic and loxodromic cases. It 
may even be viewed as a virtue of the method of [Z.l] that it adapts so effortlessly to the general case. 

It should be noted here that Guillemin was aware that the arguments of the elliptic case should extend 
to the general non-degenerate case and stated his main result, Theorem 1.4 of [G.2], for the general case. 
However, we also note that the methods used here are extensions of the methods of [Z.l], which in many 
significant respects differ from the methods of [G.2]. 

The organization of this paper is as follows: In §1 we will review the symplectic and microlocal ingre- 
dients required to construct a 'linear model' for the Laplacian near a closed geodesic 7. In §2, we introduce 
the semi-classically scaled Laplacian and the linearized Laplacian and conjugate them to the model space. 
In §3, we conjugate the resulting semi-classical model Laplacian to a semi-classical normal form to infinite 
order. In §4 we show how this semi-classical normal form induces a bona-fide quantum Birkhoff normal 
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form for the Laplacian near 7. In §5, we use the normal form to give the explicit formula (0.1) for the wave 
invariants. In §6, we show that the quantum normal form coefficients can be determined from the special 
values of (0.1) corresponding to 7 and its iterates. 

1. Preliminaries 

This section begins with a resume of the symplectic linear algebra underlying the Jacobi equation, the 
linear Poincare map, and the symplectic classification of non-degenerate quadratic forms (§1.1). It then 
summarizes the quantum aspects of the linear theory, in particular the behaviour of the quantum action 
operators (§1.2) 

§1.1: Symplectic preliminaries 



§l.la: Closed geodesies, linear Poincare maps and Jacobi fields 



Throughout this paper, 7 will denote a primitive closed geodesic of (M, g); its iterates will be denoted 
by 7" 1 . 

The space J^r of (real) orthogonal Jacobi fields along 7 is then the real symplectic vector space, of 
dimension 2n, of solutions of the Jacobi equation Y" + R(T, Y)T = (with T the unit tangent vector along 
7.) The symplectic structure is given by the Wronskian 

.(X,Y)= g{ X^Y)- 9{ §- s X, n 

The linear Poincare map P 1 is the (real) linear symplectic map on (J'J-,w) defined by PjY(t) = 
Y(t + Lj). To diagonalize it, we also complexify it as a complex symplectic map P^ on the space J^r ® C 
of complex orthogonal Jacobi fields. Here, the symplectic form is extended to the complexified space as 
a complex bilinear form oj c . Since P^ G Sp(J^~ <g> C,w), its spectrum a(P^) is stable under inverse and 
complex conjugation: thus, if p e cr(P^), then also p, p^ 1 G cr(P^). As mentioned above, we will assume 
that P^ is non- degenerate in the following strong sense: 

(l.la.l) p™ 1 ... p™» = l=> mi = (Vj,m t eN). 

In particular, the eigenvalues are simple and ±1 ^ cr(P^). 

The eigenspace of P^ of eigenvalue p will be denoted by J^ c {p) C J7^"®C We then have the symplectic 
orthogonal decomposition 

into the 'non-compact' symplectic subspace 

(l.la.3.nc) Jt C = ® P :\ P \^iJ^ C {p) 

where P^ does not belong to a compact subgroup of Sp and the 'compact' symplectic subspace 
(l.la.3.co) J^ C = ® pM=l J^ c { P ) 

where P^ does belong to a compact subgroup of Sp. This and the following decompositions are described in 
more detail in Klingenbcrg [Kl] , but also somewhat differently since in [Kl] the symplectic form is extended 
to the complexification as a sesquilincar form rather than as a complex bilinear form. 
The non-compact subspace has the further symplectic orthogonal decompositon 

(l-la.4) jt C = [J^ C {P) © J^Hp- 1 )] 

p=M<i 

into symplectic complex 2-plancs. We may rewrite this decomposition in the form 

(l.la.5) j£ C = J s ± ' C <BJt' C 

where 

7-_L,C _ ST\ 7--L.C/ _\ -T-J-.C _ n\ 7--L.C 

p:\p\<l P--\P\>1 



are the symplectically dual stable, resp. unstable Lagrangean subspaces. 
The compact subspace has the further symplectic decompositon 

(1-1-6) Jt C = J^ C {P)®J^ C {P). 

p:p—e la ,a£ (0,7r) 

Any choice of one p from a pair {p, p} determines a splitting of into a pair of dual Lagrangean subspaces. 
On the level of real symplectic spaces, we have the closely related P 7 -invariant symplectic decomposition 

(i.ia.7) = Jl e j: e Jt 2p 

into the stable, unstable and center stable real subspaces of dimensions r, r, 2p respectively. By definition, 

(l-la.8) Jl= J^{p)® J,^(p) 

peR,\p\<\ peC-M,\p\<i 

where 

P 7 |^, R(p) =p (peR) 

respectively 

_ , _,. / cosv sinv \ -i,4-iv « 

P 7 \ 7 ±,K (o) =e M . p = e fj,,v>0. 

" J -i (P) \ —sinv cosv J 

In the latter case, 1 7 7 L ' C_R (p) is the real symplectic 2-plane whose complexification equals J^-' C {p) 
J^'^p -1 ). Similarly for the case of the unstable subspace. In the center stable case, 

(l-lo-9) ^ = ©^ ±,R («i) 

where 

, / cosa sina 

"^i \ a ) y —sma cosa 

and with J^-' R {a) the symplectic two plane whose complexification equals j7 7 L,c (p) ®J^~ ,c {p) with p = e M . 

We will put r = q + 2c where <? = #{p g K, |p| < 1} and where 2c = #{p e C - R, \p\ < 1} and say 
that 7 has type (p, g, c) if it has p pairs of conjugate eigenvalues of modulus one {e M , e~ M }, g pairs of real 
inverse eigenvalues {e A ,e~ A }, and c quadruples of non-real complex eigenvalues {e ±M+±J!/ }. Here, we have 
assumed the real eigenvalues are positive for brevity of notation. Finally we may write: 

(1.10.10) jj- = j; e e j-f 

where J* = J^I is the elliptic (or real center stable) subspace, where 

(l.la.ll) J*= J^{p)@J^{p-') 

peR,\p\<i 

is the real hyperbolic subspace and where 

(l.lol2) J* = J^\p)@J^{p- 1 ) 

peC-M,|p|<i 

is the complex hyperbolic (or loxodromic) subspace. 

§l.lb: Jacobi eigenvectors and Wronskian matrix 

As mentioned in the introduction, the intertwining operator to the quantum normal form will involve a 
certain Wronskian matrix of the Jacobi equation. Roughly speaking, it is the real symplectic matrix whose 
entries are given by the real and imaginary parts of the Jacobi eigenvectors and their time derivatives relative 
to a normal frame. 

More precisely, we fix a symplectic orthonormal basis of Jacobi eigenvectors as follows: 
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complex subspace 


eigenvectors 


normalization 


elliptic plane 


P 7 Y 7 e = e ia iY?, P 7 F e j = e- ia iY e j 


uj{Y^Y^) = 1. 


real hyp. plane 


P 7 Y+ = e^Y+, P 7 Yr = e~^Yf 


w{Y?,Yr) = \. 


cx.hyp.4-plane 


P y Y+ + = e ^ +lv Y+ + , P 7 Yr- = e'^Yf-, 


= 1 




p 7 y+- = e"-'"^-, P-yYr+ = e -"+ tl, Y- + , 


lj(Y-+,Y+-) = 1 



The normalization makes sense since complex 2-planes spanned by eigenvectors corresponding to inverse 
eigenvalues are symplectic. 

We now fix a parallel normal frame e(s) := (ei(s), e n (s)) along "f\[o.L) an d denote by (Y,ej) the 
Riemannian inner product of a vector Y along 7 with the the jth normal vector. Corresponding to the 
splitting of J^j- into its elliptic, hyperbolic and and complex hyperbolic (real) subspaces, we then get a real 
symplectic 2n x 2n Wronskian matrix 

(1.16.1) W = [W e \W h \W ch ] 

formed by the 2n x 2p elliptic Wronskian matrix 

Re(Y?, ej ) Im(Yf, ej ) 
ReiY^ej) Im(Yf, ej ) 

the 2n x 2q real hyperbolic Wronskian matrix 



(1.16.2e) 



W e (s) 



*=1|...ip;j=1). 



(1.16.2/i) 



W h (s) := 



{Y^e 3 ) (Yr, ej ) 



y <Xi ,ej) <Xi ,ej) J 
and finally the complex hyperbolic 2n x 4c Wronskian matrix 

' Re(Y++, ej ) Im(Y++, ej ) Re(Yr 



(1.16.2c/i) W ch (s) := 



Re{Yi ++ ,ej) Im(Yi 



Im(Y: 



, ej) Re(Yi ,ej) Im(Yi ,ej) 



• >c;j=i,- 



To see that W is indeed symplectic, and to better understand its properties, we reconsider the Ja- 
cobi equation and Poincare map from the Riemannian viewpoint. Thus, we let V denote the Riemannian 
connection, and recall that it determines a horizontal subbundle of T(S*M) complementary to the vertical 
subbundle of the projection it : S*M — > M. Together with the symplectic structure, we get a splitting 

T(S*M) = H@V@f 

where T is the real span of 7, and H © V is the horizontal plus vertical decomposition of the kernel of the 
contact form a = £ ■ dx (or equivalently, of the symplectic orthogonal of T and the cone axis). The subspaces 
H, V are symplectically paired Lagrangean subspaces of T(T*M). Given a vector X G ^j(t)> we denote by 
X h the horizontal lift of X to -& 7 (t) and by X v the vertical lift to V^ t y The correspondence 

Y{t)^{Y{t)\Y{tf) 

then defines an isomorphism between the spaces of Jacobi fields and geodesic flow invariant vector fields 
along (7(t),7(i) (cf [Kl, Lemma 3.1.6]). That is, 

dOy Q)MQ)) : (Y(0) h ,Y(0y) - (Y(s) h ,Y( S y) 

where Y(s) is the Jacobi field with the given initial conditions. Moreover, since G* is a Hamiltonian flow 
dG s is a linear symplectic mapping from (H © t^) 7 (o), 7 (o) to (H © ^0 7 (s), 7 (s)- 

Relative to the basis {ej(s)}, the Jacobi equation is equivalent to the linear system jfc(Y, P) = JH(Y, P), 
where: (i) P = SOL, (ii) J is the standard complex structure on R 2n , and where (hi) 



H 



K 




with K the curvature matrix and with / the identity matrix. Moreover the basis {ej(s)} induces a moving 
symplectic frame {e^(s), ej(s)} of (H © ^) 7 ( s ), 7 (s)- The evolution operator for the linear system is then just 
dG s expressed as a matrix relative to the moving symplectic frame. 
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Now consider the above basis {ReY? , ImY? , Y^ ,Yj~ ,ReY^ ++ ,ImY+ + \ReY- ,ImY- } of formed 
by the eigenvectors of P 7 . By construction, it is a symplectic basis relative to the Wronskian form u). Hence 
the pairs (Y h (s), P v (s)) consisting of these eigenvectors and their time derivatives form a moving symplectic 
basis of H © V. Expressed in terms of the frame {e^(s), ej(s)} we then get a symplectic basis of M 2 ™ relative 
to the standard symplectic structure. The Wronskian matrix W(s) is just the matrix whose columns are 
formed by these basis elements, and it is therefore symplectic for each s. 

Consider now the monodromy aspect of W(s), i.e. its transformation law under time translation s — > 
s + L thru one period. Let Yi denote one of the complex eigenvectors of P 7 . Then the matrix element 
(Yi(s), ej(s)) (or with Yi in place of Yi) satisfies 

n 

(1.16.3) {Y i (s + L),e j (s + L))=p i J2tjk(Yi{s),e k {s)) 

fe=i 

where T := (tjk) is the holonomy matrix, 

n 

ej(s + L) = ^t jk e k (s). 

fe=i 

It follows that 

(1.16.4) W(s + L) = W{s)T*P- ( . 



§l.lc: Symplectic equivalence of quadratic Hamiltonians 

Let (M 2 ™, u>) be the standard symplectic vector space, endowed with linear coordinates x = (qi, . . . , q n ,pi, . . . ,p n ) 
such that oj — ^2dqi A dpi. A quadratic (real) Hamiltonian is by definition a quadratic form 

H(q,p) = l;(Ax,x) = ^lj(JAx,x) 
where (•) is the Euclidean scalar product, where A is a 2n x 2n real symmetric matrix and where J = 
Then J A G sp(M. 2n , u>), and hence its spectrum decomposes into purely imaginary pairs 



/ 

(ia, —ia), into real pairs (A, —A), and into complex quadruples (±fi±iv). We will assume, as above, that the 
eigenvalues are simple and not equal to ±1. Then H(q,p) decomposes into sums of terms of the following 
types of quadratic Hamiltonians, or classical 'actions': the elliptic type 

(l.lcle) I e {qi,Pi)~\u(,ql+p\), 
the real hyperbolic type 

(l.lclfc) I h (quPi) ■= \\qiPu 

and the complex hyperbolic (or loxodromic) type 

(lAcAch) i ch (qi,Pi,q2,P2) = ^(qiPi + qiPi) + -^>{q\vi - qiv\)- 

Note that 

q\P\ + = Re(qi + - IP2), Q2Pi - q\P2 = Im(qi + «g 2 )(pi - ip2) 

and 

{qiPi + q2V2, q\P2 - q2Pi} = 

To unify these expressions, we observe that they all have the form ^Re sa*a where seC and where a*, a 
denote symplectically dual complex linear coordinates. Indeed, in the elliptic case, a = q\ + ipi,a* = 
qi — ipi , s = a e K; in the real hyperbolic case, a = q\ , a* = p\ , s = A e R; and in the loxodromic case, 
s = n + ii>, a = (qi + iq 2 ),a* = (pi - ip 2 ). 

§1.2: Microlocal preliminaries 



As mentioned above, the wave invariants only involve the metric and Laplacian A in a tubular neigh- 
borhood of 7. In fact, as discussed in [G.l][Z.l] they only involve the microlocalization of A to the conic 
neighborhood 

(1.2.1) \y\<e, \ V \<ea 

of T*S\ — Q in T*(S 1 x M"). Here,(s, a, y, 77) denote the symplectic Fermi coordinates and ip denotes a smooth 
homogeneous cut-off function on T*(S^ x E") — which equals 1 in some conic neighborhood V of T* S\ — 
and vanishes identifically off of some slightly larger conic neighborhood. 

As in the case of elliptic closed geodesies, to put A into a microlocal (quantum Birkhoff) normal form 
around 7 ~ S 1 is first of all to conjugate it to a distinguished maximal abelian subalgebra Ay of the algebra 
iff* (S^ x E") of pseudo-differential operators on the model space S\ x E™. This distinguished subalgebra will 
depend on the type of the geodesic 7. Roughly speaking, it will consist of the tangential operator D s := i 
together with an appropriate set of quantized quadratic normal forms or 'action operators' in the transversal 
directions. 

§1.2.1: The model algebras 

To specify this 'appropriate set' of action operators, we begin by recalling that the Schrodinger repre- 
sentation of the (complexified) Hcisenberg algebra h„ <x> C on the transverse space L 2 (R n ), is generated by 
the self-adjoint operators Yj = "multiplication by yj" and by Dj — j^p- Equivalently it is generated by the 
creation/annihilation operators Yj + iDj,Yj — iDj. For our purposes, however, it will be more natural to 
choose a different set of generators depending on the (q,p,c) type of the closed geodesic. 

Corresponding to the 2p dimensional elliptic symplectic subspace we will use as generators the above 
(elliptic) annihilation/creation operators 

(1.2.1.1e) Z j :=Y j +iD yj Z}=Y j -iD Vj (j = l,...,p) 

which satisfy the commutation relations 

[Z j ,Z k ] = [Z),Zl\=Q [Z j ,Zl] = 2S ij I. 

We would like to use the real (resp. complex) hyperbolic analogues in the hyperbolic subspaces. To 
determine the analogues we note that Zj, resp. zj arc the Weyl quantizations of the symplectically dual 
(modulo a factor of 2) complex linear coordinates Zj := yj + irjj, resp. zj := yj — irjj. We use the 'dagger' 
notation rather than the adjoint notation Z* to emphasize that the dual operators are symplectically dual; 
they are also adjoints of each other, but this property will not extend to the hyperbolic cases. Indeed, 
corresponding to the 2q dimensional real hyperbolic subspace, the natural generators are the hyperbolic 
annihilation/creation operators 

(1.2.1.1ft) Yj, D yj (j=p+l,...,q) 

which of course are also symplectically dual. And corresponding to the 4c dimensional complex hyperbolic 
subspace, we the natural generators are the complex hyperbolic annihilation/creation operators 
(1.2.1cft) 

Wj := Yj+iY c+j , W] := D y -iD Vc+j , W, := Yj-iY c+j , := D Vj +iD Vc+j (j = p+q+1, . . . , c). 

We note that they satisfy the commutation relations: 

\Wj,W}] = -2, [Wj,Wj ] \ = -2 

with all other brackets zero. We will not bother to renormalize the operators to be precisely dual. 

The enveloping algebra of the Heisenberg algebra is then generated by all the above annihilation/creation 
operators, 

(1.2.1.2) £:=<Z U ... Z p , Z\, . . . , z\, Y u Y g , D yi , D Vq , W u W u . . . , W c , W c , W\,Wl . . . , W\, W\ > 

and is of course the algebra of partial differential operators on E" with polynomial coefficients. We will 
denote by £ n the subspace of polynomials of degree n in the generators. The microlocalization of this 
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algebra is the isotropic Weyl algebra W* of pseudo-differential operators on R n , in which the generators are 
assigned the order |, so that 

£ n c W n/2 

The symplectic algebra sp(n, C) is then represented in £ 2 by homogeneous quadratic polynomials in the 
generators, which have degree 1. In particular it contains the following elliptic, resp. hyperbolic, resp. 
complex hyperbolic ( loxodromic) 'action' operators: 
(1.2.1.3) 

Ij ■= Z)z h l ! ; = ^Y,l)., + D Vj Yj), lf' Re = l -Re{w}W + (WjWj)*), lf Jm = ImWjWj. 
The complex hyperbolic action operators can also be written in the form 

(1.2.1.4) 1?'** = l -{Y 3 D yj + D V] Yj + Y j+c D yj+c + D Vj+c Y j+c ), lf' Im = (Y 3 D V]+c Y 3+c D yj ) 

where the coordinates are indexed so that the dyj A dyj +c A drjj A dr/j+c-planes are the P 7 -invariant complex 
hyperbolic 4-plancs. It is then natural to introduce polar coordinates rj,<f>j on the (yj, j/j+ c )-plane so that 
the loxodromic actions operators simplify to 

(1.2.1.5) lf' Re = \{rjD rj + D rj r s ), I chJm = D e . 

We now introduce the distinguished (p, q 7 c) maximal (transverse) abelian subalgebra of W, given by 

ao 1 c\ -7- ^ je je jh jh rch,Re T ch,Re T ch,Im T ch,Im 

.Z.l.OJ X pqc . — ^l 1 ,...,l p ,l 1 ,...,l q ,l 1 ,...,l c ,...,j c >. 

Together with the tangential operator we get the (p,q,c)- maximal abelian subalgebra given by 

an -i 7\ A _^ y\ je je jh rh jch.Re jch.Re jchjm rch,Im ^ 

.Z.L.I) J\ p gc . — <. U Sl l l , i p , I 1 , . . . , l q , 1 1 ,...,J C ,I 1 ,...,J C .> . 



§1.2.2: Model eigefunctions 

An orthonormal basis of L 2 (S'| / x R n ) of joint eigenfunctions of A p . q . c is given as follows: corresponding 
to the (p,q,c)- type of P 7 we can write 

L 2 (Si x R n ) = L 2 (Sl) ® L 2 (W) ® L 2 (R q ) ® L 2 (R 2c ) 

and construct the eigenfunctions as (tensor) products of the eigenfunctions on the factors. In the elliptic 
factors, the eigenfunctions are the normalized Hermite functions j q (cf. [F], or [Z.l] for a context similar to 
the one here). In the real hyperbolic factors, the action operators are the generators of the unitary dilations 
on L 2 (R) given by 

U(6)f(x) = 6\f(9x), (9 e R+). 
Their generalized eigenfunctions are the temperate distributions 

x + 2 , x_ 2 , (aSl) 
and any f € L 2 (R, dx) has the eigenfunction expansion 

f( x )= / I+{a)x + 2+la da+ I f^(a)xZ 2+la da 
Jr Jr 

with f± a := (f,x ± 2+ta ). In the complex hyperbolic (i.e. loxodromic) factors, the actions operators are given 
by the unitary dilations in polar coordinates 

U(p)f(r,6)=pf(pr,6), ( P eR+) 

together with rotations. The joint eigenfunctions are the temperate distributions on R 2 given by 

r lt - x e m \ (t e R) 



and in a notation similar to that of the real hyperbolic case a function / £ L 2 (R 2 , rdrdO) may be expressed 
in the form 

For future reference we summarize the situation in the following table. 



Factor 


Action 


Eigcnfunction 


L 2 (Sl) 


D s 


e l 


L 2 (W) 




Hcrmitc functions 7 m , m £ N™ ; 

loiy) = iiiiy) ■= e~3M 2 

lm :=C m a\ m \..a\r*lo 


L 2 (W) 




IlJ =1 y;f f ,agR r 


L 2 (R 2c ) 


jchRe jchjm 
1 3 ' 1 .i 


n c j=1 rl j ~ L e in ^,t £ M c 



As in the introduction, we put: 



i=i i=i i=i 

and note that 
(1.2.2.2a) 

with 

1 p 1 r c 

(1.2.2.26) r femnat = -(27rfc + ^(nij + -) + ^ X j a j + E ^jtj + v i n j\)- 



2. The semi-classically scaled Laplacian 

The significance of the maximal abclian algebra A pqc will appear as soon as we semi-classically 'rescale' 
the Laplacian and conjugate the principal part, the 'linearized \/A\ to its normal form. 

Let us briefly recollect this rescaling, which proceeds exactly as in the purely elliptic case [Z.l,§2]. We 
first prepare the Laplacian by putting it in Fermi normal coordinates (s,y). It is then self-adjoint relative 
to the volume density J(s, u)\ds\\dy\ in these coordinates. To simplify, we then conjugate it to the unitarily 
equivalent (1/2-dcnsity-) Laplacian 

A 1/2 := JV^AJ-Va, 

which is self-adjoint with respect to the Lesbesgue density \dsdy\. 
We thus have: 

n 

(2.1.1) -A 1/2 = J- x l 2 d s g°°Jd s J- x l 2 + Y J~ 1/2 d w 9 ij Jd yj J- 1/2 

ij=l 

n n 

= g°°d 2 + r°d s + Y .'/'•' <>:, + E rd y* + a °- 

ij — l i—1 

Semi-classical rescaling then involves two conjugations: First, by M h = multiplication by e"nr, 
—M£AM h = -{hL)- 2 g°° + 2i{hL)- 1 g°°d s + {(tiL)-^ + A 
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and then by the semi-classical dilation Tf l f(s,y) = f(s,h ^y). The complete conjugation —Tj^M^AMhTh 
results in the semi- classically scaled Laplacian 

n n 

(2.1.2) -A, = -(hLr'g^+^hLr^^ds+zihLr^+h-^Y, 9[^^) + ''"'(E I 1^) + Ww> 

the subscript [h] indicating to dilate the coefficients of the operator in the form, /^(s, y) := / (s, hiy). 
Expanding the coefficients in Taylor series at h — 0, we obtain the semi-classical expansion 

oo 

(2.1.3) A fc ~ h ( - 2+m/2) C2- m /2 

m=0 

where £ 2 = L~ 2 , C 3 / 2 = and where 

Pi i n n 
(2-1.4). & = 2L-i[i- + -(£d 2 yj E KiMViVi)] 

We will denote the bracketed operator, the 'linearized \/A' by L. It is of order 1 in the sense of pseudo- 
differential operators (using the Weyl filtration in the transverse variables) and as will be seen below is the 
principal term in the semi-classical expansion of the square root of the rescaled Laplacian. 

(2.1. A) Appendix on metric scaling 

In addition to semi-classical scaling, we have also just introduced an independent scaling, metric scaling, 
which has to do with the behaviour of objects under dilations g — ► e 2 g of the metric. As discussed in detail 
in [Z.l], the wave invariants have well-defined weights under metric rescaling and in analysing them it is very 
convenient to rescale all objects to be weightless. For instance, as discussed in [Z.l, §1.4], an w- symplectic 
basis of Jacobi fields has weight \ and its time derivative has weight — |. To render it weightless a Jacobi 
eigenficld Y should be replaced by L~^Y, Y by L^Y etc. The resulting weightless Wronskian matrix is then 
denoted by Wl- It is essentially the weightless matrix denoted Al in [Z.l]. 

To render the coordinates (y, 77) weightless under metric rescaling, we also change variables to x = L~ 1 y 
and rewrite and the £ 2 -f ' s in terms of the ai-variables. For instance, C then takes the form: 

„ -. n n 

C = % + 2 E L ^ d ^ - E LK^s)x lXj }. 

The symplectic coordinates on the symplectic normal space T*R n to R+7 will henceforth be denoted (#,£). 

For the sake of brevity we will not draw much attention to metric scaling in the various steps to come 
in the normal form algorithm. In all cases, the role of metric scaling is identical to that in the elliptic case 
[Z.l]. 

§2.2: Conjugating A^ to the model 

We now conjugate the semi-classically scaled Laplacian A^ from L 2 (N 7 ) to the model L 2 {S\ x R") by 
means of the moving metaplectic operator ju(Wt), 

»(W L )f(s,y) := n{W L {8))f(s,y). 

The motivation for this conjugation comes from: 

(2.2.1) Proposition 

c = niWDD.niwi)- 1 

where Wl is the weightless Wronskian matrix and fi is the metaplectic representation. 
Proof: 
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First, let us ignore the scaling parameter L, i.e. let us put L = 1. The right side is then equal to 
(D s + fi(W(s))* D s /j,(W(s)). To evaluate the second term, we recall that the columns of W are Jacobi fields, 
and that Jacobi's equation is equivalent to the linear system 4j(F, P) = JH(Y, P) with P = and with 



H = 



K 
/ 



Hence, the second term is jdfj,(JH) with d\x the derived metaplectic representation. But \d^{JH) = 
l/2(Y^i=i®yi ~ Y^ij=i Kij( s )lJiyj) [F]- Re-inserting L to make all objects weightless, we get the formed 
claimed above. □ 

Thus, conjugation by Wi puts the principal term C of into the simple normal form D s . This 
suggests conjugating the full rescaled Laplacian by ^(Wl) to the 'twisted model' semi-classical Laplacian 

(2.2.2) V h = n(W* L y 1 A h n(W* L ) 
which has the asymptotic expansion 

oo 

(2.2.3) V h ~ h { - 2+ ^ ] V 2 _™ 

m—o 

with T>2 = /,£>! = 0,2?i = D s . Thus, Vh is a small perturbation of D s , and one may expect that 
perturbation theory can be used to find a good normal form for the whole of V h - 

Before doing so, we must consider which Hilbert space is the natural domain for T>h- The point is 
that the conjugation has non-trivial monodromy (§l.lb) and hence the conjugate will act on functions 
transforming correctly under the monodromy group. 

We can describe the Hilbert space in terms of quantum mapping cylinders [Z.l]. First, we consider the 
holonomy aspect, put 

(2.2.4) Cf?(R x R") := {/ G C°°(M x R") : f(s + £,,«) = ^{T)f{s, u)} 

and let TCt denote its closure with respect to the obvious inner product over [0,L) x R". Note that the 
metaplectic operator /u(T) is simply 

n{T)f(u) = f{t~ l u) 

and hence that 

C^(R x R n ) - C°°(Ar 7 ) 

where the isomorphism is simply the pull-back by the exponential map defined by the frame e(s). In other 
words, expressed in terms of Fermi coordinates relative to a normal frame, L 2 (N^) becomes the quantum 
mapping cylinder of fi(T*). Let us note however that A/j and hence all the £ 2 -i ' s are invariant under under 
/x(T), so that it will play an insignificant role for our purposes. 

On the other hand, the quantized linear Poincare map /x(P 7 ) will play an essential role. Hence we 
introduce its quantized mapping cylinder 

(2.2.5) H 7 := {/ G Lf oc (R x R") : r L f = M (P 7 )/} 
and note that 

is a unitary equivalence. Hence, the natural domain for T>h is the quantum mapping cyliner of /^(P 7 ). 

In the calculation of traces, it is simpler to work in the original model L 2 (S\ x R"). Hence we will 
also consider the conjugate of T>h under a conjugation which untwists the mapping cylinder of /x(P 7 ). That 
is, we connect P 7 to the identity by a segment of the one-parameter subgroup P 7 (s) thru I and P 1 , which 
exists by our non-degeneracy assumption on P 7 . Indeed, after diagonalizing P 7 and consulting the list of 
symplectic equivalence classes of quadratic forms, we see that 



P 7 = ex P (E Ha ^ v) ) 
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where 5 f denotes the Hamilton vector field of / and where exp o S denotes the exponential map from 
sp(n,M.) — > Sp(n,R), with sp(n,M) viewed as the Poisson algebra of quadratic functions on R 2 ™. We then 
have 

P 7 (s) = exp(sE HaX ^^) 

and quantize this subgroup as 

(2.2.6) M (P 7 ( fl) ) = «<•*-.*.(,,> = c fa E?-i "'^I- w /f - + ^f 
Conjugation by /x(P 7 ) transforms 2\ into the model semi-classically scaled Laplacian 

oo 

(2.2.7) TZ h := m(P 7 )m(Wl)A^(Wl)>(P 7 )* ~ ]T M" 2+ ^^ 2 -^ 

m— o 

with 7^2 = 7/, TZs = 0, and with TZi := TZ. All coefficients of terms in TZ are periodic in s and have weight -2 
under metric rescaling. 

3. Semi-classical normal form 

We now wish to put TZh into semi-classical normal form, in the sense of [Z.l, Lemma 2.22]. This is 
the key transitional step in putting A into microlocal normal form and is the source of the connections to 
local geometric invariants. The method is essentially the same as in the elliptic case, both in method and in 
detail. We therefore present only the first two steps in the proof and refer to [Z.l, loc.cit] for the inductive 
argument. 

As in the elliptic case, we state the result in terms of the ^-operators since the trace will later be analysed 
in this model. However, most of the proof will take place in the twisted model, where the 'linearized Laplacian' 
is D s and the equations simplify most. In the following, the notation | G means to restrict to functions in 
the kernel of TZ, that is, to elements of 7^-weight zero. In the twisted model, these are simply functions 
independent of s. In the passage from the semi-classical normal form to the microlocal (quantum Birkhoff) 
normal form, the various operators will only be applied to such weightless elements. This explains the rather 
complicated statement to follow; the result is only simple and natural when restricted to elements of weight 
zero. 

(3.1) Lemma (cf. [Z.l, Lemma 2.22] There exists an L-dependent h-pseudodifferential operator Wh = 
Wh(s, x, D x ) on L 2 (S\ x W 1 ) such that, for each s e S\, 

W h (s,x,D x ):L 2 (R n )^L 2 (R n ) 

is unitary, and such that 

oo 

WlTZ h W h ~ -h- 2 L- 2 + 2h~ 1 L- 1 TZ + h^lZ™_i_(s, D s ,x, D x ) 

1=0 

where 

(i) TZ^is, D s ,x,D x ) = TZ^TZ 2 +TZ^TZ + TZ^\ , with e C°° (Si, £r 2fe ); 

(ii) TZf_ 3 {s, D s ,x, D x )\ a = TZ?_!°(s, x, D x )\ = f^If, . . . , F p , if, . . . , 1^ I{ hR \ . . . lf Re , If hIm , lf Im )\ 
for certain polynomials f 3 of degree j+2 on R" , i. e. f 3 (If , . . . , I«, if, . . . , l£, lf Re lf Re , lf Im , . . . , lf Im ) G 

(Hi) TZ™ 2fc+1 (s, D s ,x, D x )\ a = TZ™'° 2k+1 (s,x,D x )\ a = 0; 

(iv) The terms I^,I ] -,lf h are weightless under metric scalings and all of the TZ's have weight -2. 
Proof: 

As in the elliptic case [Z.l, Lemma 2.22], the operator Wh will be constructed as the asymptotic product 

(3.2) W h := KP-rT n r=i^ ° m(P 7 ) 

13 



of weightless unitary /i-pseudodifferential operators on R™, with 
(3.3) W h k := exp(ihiQk) 

and with K*Qk £ /i*C°°(S*i) <g> £ k+2 of total order 1. The product will converge, for each s, to a unitary 
operator in ^^(M™) (we refer to [Sj] for a discussion of such asymptotic products). 
We first construct a weightless Qi(s, x, D x ) £ C°°{S 1 L ) ® E\ such that 

(3.4a) e~ ihiQ in h e ihhQ ^ \ a = [-h~ 2 L- 2 + 2h~ 1 L- 1 1Z + Tzl+...}\ 

where the dots . . . indicate higher powers in h. The operator Q i then must satisfy the commutation relation 
(3.46) {[L- 1 TZ,Q i ]+Tl i }\ o = 0. 

To solve for Qi , we conjugate back to the T> 2 ^m. , s of the twisted model by n(P-y), which transforms 1Z into 
D s . The commutation relation thus becomes 

(3.4c) {[L _1 Z) S , ^(P 7 ))*<3i/i(P 7 )] + Vi}\ = 0, 

that is, 

(3.4d) L _1 <9 s {/i(P 7 )*(2i/i(P 7 ))}| - -*{X> 4 }| 

where <9 S A is the Weyl operator whose complete symbol is the s-derivative of that of A. Since (3.4d) is 
simpler than (3.4b), we henceforth conjugate everything by ;ii(P 7 )), and relabel the operators /u(P 7 )*Q/z(P 7 ) 
by Q. The resulting D's then transform under tl like operators on the quantum mapping cylinder of /i(P 7 ). 
Our problem is thus to solve (3.4d) with an operator Q± satisfying 

r L gir£= M (P 7 )Qi M (P 7 )*. 

To solve the equation (3.4d) we rewrite it in terms of complete Weyl symbols. We will use the notation 
A(s,x,£) for the complete Weyl symbol of the operator A(s, x, D x ). Then (3.4d) becomes 

(3.5a) L~ 1 d s Qi (s, x, £) = -Wi \ a (s, x, £) 

with 

Qi(a + L,x,0 = Qi(s,P 7 (x,0)- 
We solve (3.5a) with the Weyl symbol 

Qi(s,x,£) = Qi (0, x,£)+L / -Wi \ (u, x, £)du 

Jo 

where Qi(0,x,£) is determined by the consistency condition 

(3.56) Qi (L, x, £) - Qi (0, x,g) = L [ -iVi\ (u,x,£)du 

2 2 Jo 2 

or in view of the periodicity condition in (3.5a), 

(3.5c) Qi(0,P 7 (z,0) -Qx(0,x,0 = L [ -iVi\ (u,x,Z)du. 



To solve, we use that Di \ {u, x,£) is a polynomial of degree 3 in It will be most convenient to 

express this polynomial in coordinates relative to the eigenvectors of the Poincare map. In the elliptic planes, 
we use the complex coordinates Zj — Xj + i£j and Zj = Xj — i£j (j = 1, . . . ,p) in which the action of P 7 is 
diagonal. In the real hyperbolic planes we use the real coordinates (yj,f]j) — {xj,£j), (j = p + 1, . . .p + q) in 
which the real hyperbolic part of P 7 is diagonal. Finally, in the complex hyperbolic (loxodromic) 4-spaces we 
use the coordinates Wj = Xj+ix c+ j,vij = Xj—ix c+ j,u>j = ^j—i^ c+ j,u]j = £,j+i£, c +j, (j = P+<7+l, • • -p+q+c) 
in which the complex hyperbolic part of P 7 is diagonal. 

We will denote the Weyl symbols in these coordinates by their previous expressions. We also suppress 
the subscripts by using vector notation z, z, y, r], w, w, to, u>. Thus, (3.5c) becomes 

Q i (0, e ia z, e~ la z, e x y, e~ x ri, e^w, e^w, e~^ +iv uj, e'^Q) - Qi (0, z, z, y, rj,w,w, u, Q) = 
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(3.6) = L —iT>i\ (u,z,z,y,rf,w,w,u,w)du. 

Jo 2 

We now use that T>i (u, z, z, y, r], w, w, uj 7 ui) is a polynomial of degree 3 to solve (3.5c). If we put 
(3.7a) Qi(a,z,z,y,v,w,w,w,iD):= J2 q^A^^'V* 2 ™* 1 ^ 2 ^ 1 ^ 2 

|o| + |o| + |6| + |c| + |c|<3 

and 

(3.76) Vi\ (s,z,z,y, V ,w,w,Lu,Lu)du:= ]T d ha ^(s)z a f a y b ^^w Cl ^w^u>^ 

|o| + |o| + |6| + |c| + |c|<3 

then (3.6) becomes 

|o| + |o| + |6| + |c| + |c|<3 

(3.8) . =-iL 2 d ha - abc5 z a f a y b ^w c ^ c ^^ 

|o| + |o| + |6| + |c| + |c|<3 

Under the non-degeneracy assumption on P 7 , we can solve with 

(3.9) q ha - abc -M = ~iL\l e *(»-a)°+*(ci-ei)«'+(6i-fe)A+(c a - Sa )».)-i d , io . 6c _ 

since i(a — d)a + i(c\ — c 2 )v) + (b\ — b 2 )\ + (ci — c 2 )fj) = 2nik only if a = a, b\ = b 2 ,c = c and there are no 
such (a, a, b\, b 2 , c, c) in an odd-index equation. 

Precisely as in the purely elliptic case of [Z.l], we see that Qi is a pseudodifferential operator on M ra 
with the same order, same order of vanishing, and same parity as the restriction of T>i to elements of weight 
zero. We then extend it as a pseudodifferential operator of the form 

Qi G tf^R 1 )®^ 3 

on all of 7i 7 by decreeing that it commute with s. The conjugate by /x(P 7 ) then defines a unitary operator 
W^,,i € ^hiS 1 x M") satisfying (3.4a). The corresponding twisted unitary operator with exponent Qi, i.e. 
the image of W h i under conjugation by /u(P 7 ), will be denoted W h i. 

The effect of this first conjugation is precisely as in the elliptic case: Since h^Qi is of total order 1, 

h2ad(Qi) (with ad{A)B := [B, A]) preserves the total order in^f^'*'*\ and hence W h i is an order-preserving 
automorphism of the model pseudodifferential algebra. It is moreover independent of D s and has an odd 
polynomial Weyl symbol, so that 

(3.10) htad{Qx) : h%¥(M.)(g>£™ -» [tf'-^R) ® £™+ 3 + *'(R) ® £™ +1 }. 

Finally, the d±. n 's have weight -2, the variables z have weight and hence the gi. m n 's have weight 0. 
Consider now the element 

^1 ■■= W* hi V h W h i G *2(Rl x R n } 

which can be expanded in the semi-classical series 

CO .j 

(3.11) V\-Y. h ~ 2 ^ E l J^dQ h yV 2 ^ 

n—o j-\-m—n 

oo 

:= h- 2 L- 2 + h- x L- x D s + ]T ^- 2+ ^P|_ t . 

n— 3 

An obvious induction as in the elliptic case [loc.cit.] gives that 

ad(Qiyv 2 _m. G C°°(R,£™ +j - 4 )D 2 s + C°°(R,£™ +j - 2 )D s + C°°(R, £ e m+J ). 

1 

It follows that D 9 2 _n has the same filtered structure as V 2 -^. 

2 2 
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We carry this procedure out one more step before referring to [Z.l] for the inductive argument, since 
the even steps behave differently from the odd ones. We thus seek an element Qi(s,x,D x ) e ^(S 1 x R") 
and an clement f a (If, jf , ... , J* lf' Re , lf' Im , . . . , lf< Re , I^ 1 ™) g A so that 

(3.12) V\ ~Wl{D?W hl =h- 2 L- 2 + h- 1 L- 1 D s + h- 1 *v\ +Vl{s,D s ,x,D x ) + ... 

2 

with 

(3.13a) Vl(s, D s ,x, D x )\ a = f (If, . . . , I c c h ' Im ) 

with Whi — e 1 ^ 1 , and where the dots signify terms of higher order in h. Note that V\ —T>j, so that (3.12) 

2 2 

implies that 

(3.136) {hH\ + V\}\ = /„(/«, . . . ,/f < /m ). 

2 

The condition on Qi is then 

(3.14a) {[£> a , Qi] + Vl}\ = f {It, I c c h ' Im ) 

or equivalently 

(3.146). dsQ^o = {-VI + f (If, . . .,I c c h < Im )}\ 

We solve (3.14b) by again expressing everything in terms of complete Weyl symbols relative to the 
eigenvector coordinates. Thus we rewrite (3.14b)) in the form 

L~ 1 d s Qi(s 1 z,z, y,r), w,lo,w,lu) = 

(3.15a) = -i{Vg\ (s, z, z,y,rj,w,w,w,u>) - / (|zi| 2 , ■ ■ ■ , \z p \ 2 , yirji, . . . ,y r Vr, Rewuj, Imwuj)} 

or equivalently 

Qi(s,z,z,y,r],w,w,w,u)) = Qi(0, z,z,y,r],w,w,w,u)) 

(3.156) -iL / [Dg \ (u, z,z, y,rj, w, w,w,u>) - fo(\zi\ 2 , ■ • • ,\z p \ 2 ,yirji, . . . ,y r ri r ,Rewuj,Imwuj)}du 
Jo 

and solve simeltaneously for Qi and f a . The consistency condition determining a unique solution is that 

Qi(L,z,z,y,r],w,w,w,u) = Qi(0, z,z,y,r],w,w,w,u) 

(3.16a) -iL / [D£\ (u,z,z,y,r),w,oj,w,Lj) - f (\zi\ 2 , . . . ,\z p \ 2 ,yir]i, . . . ,y q ri q ,Rewuj,Imwoj)}du. 
Jo 

or in view of the twisted periodicity condition 

Qi (0, e ia z, e- la z, e x y, e~ A ?7, e^w, e^ iv uj, e'^w, e'^Q) - Q x (0, z, z, y, rj, w, oj, w, Q) = 

(3.166) ~iL{ T>g\ (u,z,z,y,T],w,ui,w,uj)du- Lf (\zj\ ,y j r)j : Rew m uj mi Imw n uj n )}. 

Jo 

i 

In the spirit of the previous step, we use that T>£ \ a (u, z, z, y, 77, w, u>, w, u>) is a polynomial of degree 4 to solve 
the equation. We put 

gi(s,z,z,y,77,w,w,w),w) = q U aabcd s )z a z a y bl ri b2 w Cl LU C2 w El u> 52 , 

|o| + |o| + |6| + |c| + |c|<4 

and in an abbreviated notation, 

(3.17) f (\z\ 2 ,yr],Rewu,Imwuj) = ^ c oMn \z\ 2k (y ■ V ) l {Rew ■ u) n ^{Imw ■ cu) n \ 

|fe| + K| + |n|<2 

and 

vl\ (s,z,z,y, V )du:= ]T dl a - abc£ ( S )z a z a y b ^ b2 w c ^ C2 w c ^ E2 , 

|o| + |o| + |6| + |c| + |c|<4 
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and finally 

-i 1 f L i 

^l-.aabcc : = Yj J ^o:aabcc( S )^ S 

As above, we can solve for the off-diagonal coefficients where either a ^ b or m ^ n 

(3.18a) gi;aB5cc(0) = -lL 2 (l - e i(a-a)a + i(c 1 -c 1 )-+(6 1 -6 2 )A)+(c 2 -c- 2 ), t) -ljl a _ bc _ 

and must set the diagonal coefficients with a = a, c\ = c\,b\ — b 2 , c 2 — c 2 equal to zero. The expression in 
(3.18a) is well-defined by the non-degeneracy assumption. The coefficients c Q ke are then determined by 

_i 

(3.186) c oUn = dl kkUnn . 

It is evident that Q\ and fo(If, ■ ■ ■ ,Ic h ' Im ) are even polynomial pseudodifferential operators of degree 
4 in the variables (x, D x ), that Qi is weightless under metric rescalings and that the coefficients c ke n are of 
weight -2. 

The rest proceeds as in the elliptic case. □ 

4. Normal form of the Laplacian: Proof of Theorem I 

We now use the semi-classical normal forms to put the Laplacian into quantum BirkhofF normal form. 
Essentially this amounts to taking direct sums (or integrals) of the semi-classical normal form over various 
internal Planck constants. 

Proof of Theorem I: As in the elliptic case, we make the transition from the semi-classical normal form 
to the quantum Birkhoff normal form by using genreralized eigenfunction expansions for the model algebra. 

From the table in §1.2.1 we see that a function / E L 2 (Sl x 1? x R« x M 2c e ) can be expanding in terms 
of joint „4j, gc -eigenfunctions as: 

f(s,x,y,r,6) = J2 E / / f±(k,m,n,a,t)e ir ^e i<n ' e > lm (x)y^-K it - 1 dadt. 

± (K,m,n)GN 1+p+e 

Here as in §3, we have used the notation x for linear coordinates on the elliptic factors, y for those on the real 
hyperbolic factors, and polar coordinates Wj — rje l9j in each -plane of the complex hyperbolic factors. 
We also employ a multi-index notation. 

We now assemble the semi-classical intertwining operators into the Fourier-Hermite-Mellin -series- 
integral intertwining operator 

(4.1) W 1 : L 2 (Sl x R p x W x K 2c , dsdxdydw), ^> L 2 (Sl x R p x R q x M 2c , dsdxdydw) 



W -<Y. E / / f(k,m,n,a,t)e ir ^ s e l<n > e> lm (x)y± K lt - 1 dadt = 

± (k,m,n)eW+P+' jRQjm+C 

EE / / Hk^^^a^^-^Wk^rte^^jMyrK^dadt 



with 



w kmnrt ■= n{w{s)*)w r -, news)* 

kqnat 



Also, the dilation operators will be assembled into the operator 
(4.2) T : L 2 {Sl x W x R q x R 2c , dsdxdydw) -» L 2 (Si x R p x R q x R 2c , dsdxdydw)), 

T E E II f(k,m,n,a,ty r ^ s e l<n - e> lm (x)y±- h r lt - 1 dadt = 



± (fc,m,n)eN 1 +J > 



= E E II />.™.»>M)r$!l +|t|) e^^ 

± (K.m,n)(!N l +P+ c 

Here, we used that the hyperbolic eigenfunctions are eigenfunctions of dilation operators. 
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It follows, formally, from the semi-classical normal form and from the eigenfunction expansion that 



(4.3) W~ 1 T- 1 kTW 1 ~£ 2 + f a (If, ...af/ m ) + 



/i(Jf,-,j£- Jm ) 



We now show that the intertwining operator is actually a standard Fourier Integral operator (in the Weyl 
operator, or isotropic, sense) and that (4.3) holds modulo the the kind of error stated in Theorem I. 

The proof is again similar to the elliptic case, so we concentrate on the novel aspects and refer the 
reader to [Z.l, Proposition 3.4] for the remaining details. As before, we will not be as careful here as in [Z.l] 
to express things in weightless terms relative to metric rescalings. 

(4.4) Propostion TW^T^ 1 is a (standard) Fourier integral operator, well-defined and invertible on the 
microlocal neighborhood (0.1) in T* (S\ x R"). 

Sketch of Proof: 

We first consider the unitarily equivalent operator TWT' 1 in the microlocal neighborhood (1.2.1) in 
the twisted model, with 

(4.5) W : H a - H a 

W(e irk ™™ tS e i<n > 6> lm (x)y±~ V*" 1 ) := e ir " m ™ tS W-i e l<n ' e> lm {x)y± "M" 1 , 

kmnat 

and with f the dilation operator analogous to (4.2) but relative to the basis e lTkmnatS e l<n,e> j q (x)y± 2 r '* _1 . 
We then factor fwf- 1 as the product TWT 1 - 1 = ffvf' 1 where: 

(4.6) V : H a -> Lf oc (R x R x E") 

V = U^ ex P [iD; i Q i (s',y,D y )] 

that is, 

Ve ir ^'e i<n ' e> 'r m (x)yt~ i r it - 1 ) := e ir '»-¥ r -, (s', x, y, w, D Xl D Vl D^e^^jMy^'K^ 1 , 

kmnat 

and where 

(4.7) j* : C°°(M x K x E") — > C°°(R x K") 

j*.f{s,x) = f(s,s,x) 

is the pullback under the partial diagonal embedding. 

The discussion of V and the proof that j*V is a standard Fourier Integral operators goes precisely 
as in [Z.l, Proposition 3.4]. The effect of the dilation is to convert the isotropic calculus into the pure 
polyhomogenous calculus (see also [G.l] for this aspect) and then the power of D s insures that the phases 
are all homogeneous of degree 1 and vanishing to higher and higher order along 7 (by one step as the index 
j increases by one unit). Hence the phase of the infinite product has only finitely many terms of a given 
vanishing order and converges as a formal power series in the transverse variable. A convergent product can 
be defined (by Borel summation) of the phase (cf. [Sj]). 

The proposition then follows by expressing 

TW-yT- 1 =T / j,(Wyf- 1 TWf- 1 fn(W)T- 1 

and noting that T/i(>V)T _1 is also a standard Fourier Integral operator. □ 
We now complete the proof of the quantum normal form Theorem I for \/A, stated in an equivalent 
form in terms of W 7 . As in the introduction, the notation A = B means that the complete (Weyl) symbol of 
A — B vanishes to infinite order at 7 and Oj^ m denotes the pseudodifferential operators of order m whose 
Weyl symbols vanish to order j at (y,rj) = (0,0). Here, pseudodifferential operator can refer to cither the 
standard polyhomogcneous kind, or to the mixed polyhomogcneous- isotropic kind as in ^(S^) (g> W , in 
which case the total order is defined to be m = k + l. To simplify notation, we will denote the space of mixed 
operators of order m by ^mxi^L x R"). 
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(4.8) Lemma Let TW 7 T 1 be the Fourier Integral operator of Proposition (4-4), defined over a conic 
neighborhood of R + j in T* (S^ x R"). Then: 

W-'T-'VATW^ ee P^C, If, ...,lt Im ) + Po(C, If, ...,/ 2 c c Mro ) + - • • mod ®^0 2{m+1 _ k) ^- x k {Sl xR"), 
where 

, r n^(T e jchjm-, JSJ/re jchjm^ 
(a o^ TD ( f T ch,Im\ —fx VI \ 1 1t--i 1 2c ) , 1 J 2 K 1 !' J 2c /> , 
l 4 - y J nl*-! 1 1 > •••) -'2c J = Z - H ^£ 1 (TC) 2 

oo [/c— m] / t C T ch,Im\ 
P /r re T ch,Im\ _ Pfc V-*!; •■•;- / 2c J 

fc=m ^ ' 

wit/i p^ 5 m ' ; /or m= -1,0,1,..., homogenous of degree l-m in the variables (If, ...,L c ^ lm ) and of weight -1. 
Proof: 

As a semi-classical expansion in the "parameter" h = j^, (4.3) may be rewritten in the form : 

(4.io) w~ 1 t~ 1 Vat w 7 ~ £ + pi(AW 2 c c Mm ) + vM^0^1 +.... 

From the fact that the numerators fj{If, ...,I^ ,Im ) in (4.3) are polynomials of degree j+2 and of weight -2, 
the numerators Pk(If, ■•■7^2c' / " 1 ) are easily seen to be polynomials of degree k + 1 and of weight -1. Hence 
they may be expanded in homogeneous terms 

/a -i -i \ [fc+1] , [k] . [o] 

(4-11) Pk=P k +Pk +---Pk- 

with p£ the term of degree j and still of weight -1. The right side of (4.12) can then be expressed as a sum 
of homogeneous operators: 

(4.12) Pi (£,/!, ...,i£ lm ) + P {C,I{, ...,i£ lm ) + ... 

with 



J 2 ] (Te T ch,Im\ J3]/7-e jchjm, 

(aiq\ r> I r re T ch,Im\ _ f . Pi \ 1 1t---i 1 2c ) , V2 \ 1 1t--i 1 2c I 

(4-13) l J i(L,l 1 ,...,l 2c ) = L-\ — 1 (LCf 

oo [k-m],j e jch,Im\ 
P I C T e Tch,Im\ _ Pk \ 1 l'---' I 2c ) 

m\ i 1 ' 2c J - (L£) fc 
fe=m 

We claim that: 

(re jch.Im\ / j e jch,Im\ / j e T ch,Im\ 

(4.14) w~ 1 t~ 1 Vat w 7 - [£ + } + M ^ffl } + • • • + Pm[ \ L ' c ^ c ] ] 

e ©rJo 1 o 2(m+ i- fe )^(^i x R"). 

Indeed, from the analysis of the remainder terms in the semi-classical normal form (see Lemma (3.1 (i)) 
and [Z.l, Lemma 2.22]), we have 



(4.15) Pi(C,It, ...,lf/ m ) ~[C+ ^ i + • • • + 



J2]/7-e jch,Im, \N+l], Je jch.Im, 



L£ (L£)* 
G 02(^+2)^(51 x R") 



and also 



N J. k - m \(Te jch,Im\ 

(4.16) P_ m (£, If, I 2 t' /m ) - £ ( S 2C e 2(J v + i- m) * m ?(5l x R"). 

k=m 

Hence the expansion (4.10) is also asymptotic in the sense of = . For the statement of Theorem I in the 
introduction, it is only necessary to conjugate under niW). The rest proceeds as in the elliptic case. □ 
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5. Wave invariants and residue trace: Proof of Theorem B 



The purpose of this section is to show that the wave invariants have precisely the same relation to the 
coefficients of the quantum normal form in the non-degenerate case that they have in the elliptic one. The 
characterization of the wave invariants in Theorem I will then follow from Theorem A of [Z.l]. 

We will need to use some further notation and results from [Z.l]: First, the kth wave invariant of a 
positive elliptic operator P at a non-degenerate closed bicharacteristic 7 will be denoted T lk (P) . According 
to [Z.l, Proposition 4.2] we then have: 

(5.1) MP) = r 7fe (Pp fc+4 + P^ 2 + • • • + P° fe _i) 

where Pp k denotes the first k terms in the Taylor expansion of the jth homogeneous part of the complete 
symbol of P at 7. Thus, r 7 fc(P) involves the (2k+4)th jet of the principal symbol, the (2k+2)-jet of the 
subprincipal term, up to the zero-jet of term of homogeneity order (-k-1). 

As in [Z.l, §4], we will also rewrite the normal form in terms of of D s and H^ x.^.u) using that 

[LTIY {LD s y K LD S 2 { ){ LD S ' >' 

By (5.1), we can drop the DJ^ k+Vl+v H^^i^v) anc ^ higher terms, so T>k+\ can be written in the form 

(5.2) V k+1 = LD S + H aM ^ + — + ^—y 2 + • • • + {LDs)k+1 

modulo terms which make no contribution to T kl . 
We then use the fact ([Z.l, (4.3)], [Z. 2]) that 

(5.3) r lk {P) = resD k t MDs,y,D y )e ltp \t=L 

where res is the non-commutative residue and where ip t (D s ,y,Dy) denotes a microlocal cut-off to the cone 
(1.2.1). Note that in contrast to the elliptic case, the microlocal cut-off cannot be constructed in A Pt q tC since 
the neighboorhoods given by 7 < ea in terms of mixed hyperbolic-elliptic actions are of infinite transverse 
symplectic volume. This does not pose a genuine problem, but accounts for a number of modifications to 
the elliptic case in [Z.l]. For the gauging elliptic operator we use LD S . 
To simplify the notation we will put 

~ l re jch,Im\ ~ ( re T ch.Im\ ~ / Te T ch.Im\ 

(5.4) r k+1 (D 8 ,I l7 ...,I 2c ).= — + ^—j- 2 + ... + jzd^+i • 

so that: 

(5.5) r 7fc (VA) = P e ^ = oTr79^ e (^,y,^)e 4 * [ * (2 " iD3+ff - A ' ( -^ )+Pfc+1 k^)" z |*=L- 
As in the elliptic case, a key role will be played by the (formal) trace 

T(a, A, (fi, v)) := Tre iH <*^»> . 

Its precise definition is the following: Since e lHa - x -^^ is an element of the metaplectic representation \i of 
Mp(n, K), T(a, A, (fi, v)) may be identified with the character Ch of p, evaluated at the associated element 
P 7 = exp(E.H a x (,,„)) ^ Mp(n,R). Here, H a _x.(^.u) denotes the quadratic function on M. 2n which gives the 
complete Weyl symbol of the corresponding action operator, and as above exp o S denotes the flow at time 
1 of its Hamilton vector field, or, more correctly, the lift to Mp(n,M.) which corresponds to e rH <*^,(»,v) under 
(i. 

Since Mp(n,M.) is a semi-simple Lie group, the character Ch is a real analytic function on the open 
dense subset Mp(n, R) re9 of regular elements of Mp(n, R), where it is given by the Harish-Chandra formula 
[Kn]. We will need below the explicit formula for Ch{x) in terms of the eigenvalues of x. For elements of 
Mp(n, R) not having 1 as an eigenvalue, we recall that Ch(x) is given by 

Ch(x) = 



y/\det(I - x)\ 
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where a is a certain Maslov index. For non-degenerate x with p pairs of eigenvalues e ±laj of modulus one, q 
pairs of positive real eigenvalues e ±Xj and c quadruplets of eigenvalues e ± (vj±™j) j Ch{x) is therefore given 
(up to a Maslov factor) by 

(5.6) T(a, A, (p, i/)) - Ig. - • n?_j j- • — . 

Here we have selected one eigenvalue p from each symplectic pair p,p~ x (see §1.1-2). The ambiguity is fixed 
by the Maslov factor i' 7 , which can (and will) be ignored below for the sake of brevity. 
We can now give: 

Proof of Theorem B: Since e 2mLDs = I on L 2 (S l L ) we have 

afe 7 = T 7 fc(\/A) = 

(5.7) Re Sz=0 TrMD s ,y,D v )[j2(TTLD s + H a>XM ) + V k+1 ] k e iH ^^ e iLV ^ (LD s y z . 

In view of the microfocal cutoff, the operator under the trace is of trace class for Rez sufficiently large. Indeed, 
in estimating the trace we may eliminate the unitary factors and we are then left with a pseudodifferential 
operator whose complete symbol is a polynomial in (a,y,r]) times a factor of <j~ Rez ip e (a,y,ri). The integral 
in the transverse (y, r]) variables is bounded by the volume of the ball {y 2 +i] 2 ) < a and hence is of order a n . 
Since a pseudodifferential operator is Hilbert-Schmidt if its Weyl symbol is in L 2 , the operator under the trace 
is Hilbert-Schmidt for Rez > n and in particular is of trace class. Moreover, since it is the non-commutative 
residue of a Fourier Integral operator, one knows apriori that it admits a meromorphic continuation to C 
with at most simple poles [Z.2]. Hence the residue is well defined. 

As in [Z.l], we view the trace as a function of the parameters (a, A, /i, v) and use the explicit form of 
the exponential in e lHa - > -' ( ^' ,y > to rewrite (5.7) in the form 

oo V Q 

Res z=0 J2 n- z TrMDs,y, D y ){[-(2nn + a 3 D ^ + A A 3 + 

n=l *~ 3=1 j=l 

2c 

(5.8) +X>;0 W + VjD Vi ) + n+i(n, D ai ,..., D V2c , L)] k e iLVh+1 ^ n ' Dai ,-^ 2c ,D e ^, M ^>}. 

Here, we have used that D s commutes with (y, D y ) to replace it by its eigenvalue in the s-trace, and we have 
repeatedly used identities of the form F(D x )e lxP = F(P)e lxP (x e R). 

Since Vk+i(n, D ail . . . ,D V2c ,L) is a symbol of order — 1 in n with coefficients given by polynomials 
in the operators D aj (etc.), we can expand the kth power in (5.8) as an operator- valued polyhomogeneous 
function of n. At least formally, we can also expand the exponential e lL ' p k+i(n,D ai ,...,D l/2c ,L) m a p 0wer ser i e s 
and then expand each term in the power series as a polynomial in rT 1 . Collecting powers of n, the right 
side of (5.8) may be written in the form 

oo oo 

(5.9) Res z=0 J2 J2 n- z+k - j F k , k - 3 {D ai , . . . , D„ c+i „ c , )Tr^(n, y, D v )e iH "-^^ , 

n— 1 j—o 

with !Fh,k-j(D a , . . . , D IJtc+ i Vc , D ^-i^L) the coefficient of n k ~ 3 in (5.8). The expansion of the exponential is 
justified as in the elliptic case: as in [Z.l, (4.29)] we may write 

e iLVk + l := e N {iLV k+1 ) + {iLV k+1 ) N+1 b N {iLV k+1 ) 

with e^iix) = 1 + ix + • • • + ^ , with V k +\ short for V k +i{n, D ai , . . . , D V2c ) and with b^{ix) a bounded 
function. The ew term contributes a finite number of terms of the desired form (5.9). For the remainder, we 
expand (iLV k +i) N+1 as a polynomial in rT 1 with coefficients given by operators QN P (D ai , . . . ,D V2c ) and 
observe that each term has a factor of n _JV_1 . For each such term, we remove the coefficient operator Qjyp 
from the sum ^2 n , as above, leaving only the factor of &at. Since bjqiix) is a bounded function, it follows 
that bN(iLV k +i) is a bounded operator on L 2 ; and since each term of the resulting sum has at least the 
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factor n z N 1+k (possibly multiplied by a further negative power of n), we see that the remainder is a sum 
of terms of the form 

(5.10) Res z=0 Q Np {D ai , . . . , D V2c ) ]T n- z+k - N - 1 - l b N (tV k+1 (n, D ai ,... ))Tti> e (n, y, D y )e iH ^^ . 

kn 

We then observe that the sum is bounded by J2m=i m~ Rez ~ N ~ 1+k+n , hence converges absolutely and 
uniformly for Rez > — N + k + n. It follows that for N > (n + k) the sum in (5.10) defines a holomorphic 
function of z in a half-plane containing z = and since the operations of taking the residue in z and 
derivatives in a commute, each term (5.10) is zero. This justifies (5.10) and shows that it is actually a finite 
sum in j, say j < M (in fact M=(k+l)(n+k+l)). 

The residue in (5.10) is therefore well-defined and independent of e. Since Trip e (n, y, D y )e lHa ' x '<-'*-' y > — > 
T(a, A, jj,, v) in the sense of distributions as e — > oo we must have 

oo M 

(5.11) a 7fe = i?es z=0 ^2 ^ m ~ z+k ~ 3 Fk,k-]{D a , . . . , D^ c+Wc , D^-^J T(a,X,fi,v) 

m— 1 j — o 

M 

= Res z = ^2 C( z + 3 - k)Tk,k-j{p a , . . . ,D Pc+tl/c ,D^ w J T(a,X,fj,,u). 

Here, £ is the Riemann zeta-function, which has only a simple pole at s = 1 with reside equal to one. It 
follows that the only term contributing to (5.11) is that with j = k + 1 and hence we have 

(5.12) a 7fe = Fk,-i(D a , . . . , D IXc+iVc ,D IXc _ iVc )T{a, A, jj,, v). 

It follows that the wave invariants consist of the geometric data contained in the coefficients of the Tk-x^-, 
and hence in the normal form coefficients. But the algorithm for constructing the normal form is essentially 
the same as in the elliptic case and so the geometric data entering into the normal form coefficients is of 
precisely the same kind. □ 

6. Inverse Problems: Proofs of Theorem II and Corollary II. 1 

Our first goal in this section is to prove that the wave invariants of 7,7 2 , . . . determine the quantum 
normal form coefficients at 7. 

Proof of Theorem II: 

We recall that the quantum normal coefficients are the coefficients of the action monomials in the action 
polynomials Pj(If, . . . ,I% C ' m ) of (5.2). These coefficients determine, and are determined by, the coefficients 
of the monomials in the polynomials p v (If , . . . , I^' 1 ™) in Theorem B. They also corresponding bi-uniquely to 
the coefficients of the classical action monomials in the complete symbols of cither set of action polynomials. 

We also observe that the quantum normal form coefficients determine, and are determined by, the 
coefficients of the constant coefficient partial differential operator (PDO) 

(6.1fc) n,-i(D a ,D x ,D^ +tl ,,D^ lu ):= ]T C k ., abciC2 D a a D b x D c ; +lv D c ;_ lly 

(o,6,ci,C2)eN»:|o| + |6| + |ci| + |c 2 |<fe+l 

where a £ N p ,6 G N 9 ,ci,c 2 G N c . This can be proved easily by induction on k: In the case k=l, Tk-\ is 
obtained from p\ by substituting the variables D ail etc. in for the variables /f, etc. Assuming inductively 
that we have determined the coefficients of pi,..-,Pk from those of T\ . . . , T\.k, we note that pk+i 
contributes to the residue (5.7) for the first time at the k + 1st stage. Since it comes with the denominator 
D k+1 , it only contributes to the residue when composed with D k . From the form of (5.7) it is clear that 
only one term involving Pk+i contributes non-trivially, and that is the one which appears in the linear term 
in the expansion of the exponential. Hence, its contribution to Fk+1,-1 is again just the substitution of the 
variables D ai , etc. in for the variables If, etc. 

We next observe that the quantum normal form of A at any iterate 7^ of 7 is the same as for the 
primitive 7 itself. Hence the PDO Tk.-i is independent of the number N of iterations. On the other hand, 
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under the iteration 7 — > 7^, the Poincare map transforms by P 7 « — > . Therefore the expression in (5.12) 
for the kth wave invariant of 7^ is given by: 

{6.2k, N) - n ?=i (1 _ w j} ■ n *=i^T7i) ' ^ (i _ J^W) (i _ ^fW) lK>*»')=iv(a,A,^)- 

It therefore suffices to prove that for all k the coefficients of the PDO Tk,-\ can be determined from its 

values (6.2 k, N) on T(Na, NX, Nfi, Nv) for N = ±1, ±2, 

We begin the proof by noting that (6.2 k, N) can be rewritten as: 

(6.3ft, JV) rT| =1 e* jAr ^ • ITj^ei^ • n^ei^w+^e*^-*^- 

■^ fc ,_i(i>a' + + i^'w + i^'-^' + J)[nf=i- 



TT<? |_ _ TTC \ \ II 

(1 — e 3) (1 — e vf a 3 ; ) (1 — e vf a 3') 

Making the substitutions D a — » Z? Q + ^ (etc.) in (6.1 k) we obtain a new PDO whose coefficients C' k . abciC2 
correspond in a bi-unique way with the original Ck-,abc 1 c 2 ' s - Hence it will suffice to show that we can determine 
the Cfc ;a h ClC2 's from the values (6.3 k,N). 

To do so, we will regard (6.3 k, N) as the values at integral points z — N of a function of z. From the 
fact that 

D a (l - e 4 ")- 1 = [(1 - e M )- 2 - (1 - e 101 )- 1 } 

we see that this function is a polynomial in (1 - e^)" 1 , (1 - e* A ')~\ (1 - e z ^ +iv ))- 1 , (1 - e^"^))" 1 . We 
clear the denominators to obtain the entire function 

(6-4M) E ^;a6c 1 c 2 [n^ =1 (l - e" *) ■ IJ* =1 - 

(o,6,ci,c 2 )eN n :|o| + |6| + |ci| + |c 2 |<fc+l 

•(1 - e zX3 ) ■ II^_ 1 (1 - e z ^ i+iv ^)(l — e z (^^ il/ 3))]( fc + 1 ) . Il p -_ 1 ei izaj ■ lT 3 _ 1 e^ zAj ■ Yi c ._ ie ^ z ^ 3+iu ^ e^ z ^ 3 ^ iv ^ ■ 

r^a j~.b t-)Ci 7->c 2 rjrP j; tt<7 |_ ttc 1 1 11 

(1-e ') (1 — e 3) (1 — e^ p 3 j') (1 — e vf J j ; ) 

which is an exponential polynomial of the form 



(6.5ft) £ 



z03+(i,...,l),(a,A^)) 



C/c;,3e 



(6.6) Lemma 1 TTie coefficients Ck-.abc 1 c 2 can be determined from the coefficients c^p in (6.5 k). 

Proof: We first show that the coefficients Ck ia bc lC2 w hh \a\ + \b\ + \ci \ + \c 2 \ = fc + 1 can be determined from 
the Ck-ff's. Indeed, before multiplying by 

[ n j=i(! ~ e " a3 ) ■ n^iC 1 - eZAj ) ' n j=i(l - e z( ^+^ ) )(l - e *(^-^))](fc+i). 

.U p _ 1 e^ izaj ■II^_ 1 e' zAj ■ n^_ 1 e' z ^ +il/ ^e' z ('* j- " / ^ 
C^.. abciC2 is uniquely determined as the coefficient of the monomial 



jtP _ gizaj-j-faj+l) _ tt9 1 . (1 _ e zA i ) -( b j + 1 ) . IIj-i(l — ) ) -( c ij + 1 ) ( 1 _ e *(W -""4 ) )-( c 2j+l) _ 

Since we are multiplying thru by a quantity independent of a,b,c\,C2, it follows that C' k . abciC2 is uniquely 
determined as the coefficient the of the monomial 



II 3 ?_ 1 (1 - e izOl 0}( k + 1 )-( a 3+ l ) e h iza 3 . II 9 _ 1 (1 - e zA 3)( fe + 1 )-( b 3+ 1 ) e 2 



.n c _ 1 (l — e z (A l J+ ily j))(' c + 1 )-( c ii + 1 )(l _ e z (H- iv j)}( k + 1 )-( c 2J+ 1 ) e h z (H+ i ' / j) e h z (H— i ''j). 
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Expanding into an exponential polynomial, we find that C' k . abciC2 is uniquely determined as the coefficient 
of the monomial 

n^_ ie w((fc+i)-(oj+i)+3)«j . n]_ 1 e z((fc+1)_(bj+1)+ i)^ • 
.n^_ 1 e z ((' i; + 1 )-( c i3+ 1 )+^)(w+ i ^)e z(( ' c+1)_(c2 J )+ 3)(w- i ^)). 

Uniqueness follows from the fact that the vector (3 + \ is a minimal element of the set of exponent vectors 
occuring in (6.5 k). Since C' k . abciC2 — Ck- a bc lC2 when \a\ + \b\ + \c\\ + |c 2 | = k + 1, we have determined 

We then remove the C k;abciC2 (D a + \) a {D x + \) b (D^ +iv + \) Cl {D^ lu + \f terms with \a\ + \b\ + 
\ci\ + | C2 1 = k + 1 in (6.5 k). This leaves only terms with coefficients Ck-, a bc lC2 with \a\ + \b\ + \c\\ + \c2\ < k. 
Hence we can continue the process of recovering coefficients until the end. □ 



Let us now rewrite 
(6.6a) ^2 Ck-pe 

in the form 



2</3+(±,...,i),(a,A,^)} 



M 

ZUJj 



(6.66) X! a -J fee 



(6.7) Lemma The complex exponents ujj in (6.6b), together with tt, are independent over the rationals. 
Moreover, the coefficients Ck-p can be determined from the coefficients ajk 

Proof: The ujfs are rational linear combinations of the exponents ctj, \j, fj,j,Vj, which by assumption 
are independent, with tt, over the rationals. This independence also implies that the exponents (/3 + 
(i, . . . , i), (a, X, fi, v)) are all distinct. Hence the coefficients in (6.6a)-(6.6b) are the same. □ 

The proof of Theorem II is thus reduced to the following general statement about exponential polyno- 
mials. 

(6.8) Lemma Suppose that the exponents ujj of an exponential polynomial (6.6b) are independent (with 
tt) over the rationals. Then the coefficients ajk of can be determined from the values of this polynomial at 
z = N e N. 

Proof: If not, there would exist a polynomial with the given complex frequencies which vanished at all 
integers z — N. But the different terms e ZLOj , e zu)k have different exponential growth rates along z — N or 
z = —N (N G N) unless Reu)j = Reu>k- Let us write the large sum as a sum of smaller sums with a common 
Reu>. Each of the smaller sums must separately vanish for z e N. Multiply each one by the relevant factor 
of e~ Reul . Each then turns into an exponential polynomial with imaginary exponents, which vanishes for 
all z — N. Since the frequencies are independent (with tt) over Q, each of these polynomials must vanish 
identically if it vanishes at integral points. Hence the coefficients ajk are uniquely determined by values of 
the large sum at integral points. □ 
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